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Abstract 

We offer several perspectives on the behavior at infinity of solutions of discrete Schrodinger 
equations. First we study pairs of discrete Schrodinger equations whose potential functions 
differ by a quantity that can be considered small in a suitable sense as the index n — > oo. With 
simple assumptions on the growth rate of the solutions of the original system, we show that 
the perturbed system has a fundamental set of solutions with the same behavior at infinity, 
employing a variation-of-constants scheme to produce a convergent iteration for the solutions 
of the second equation in terms of those of the original one. We use the relations between the 
solution sets to derive exponential dichotomy of solutions and elucidate the structure of transfer 
matrices. 

Later, we present a sharp discrete analogue of the Liouville- Green (WKB) transformation, 
making it possible to derive exponential behavior at infinity of a single difference equation, 
by explicitly constructing a comparison equation to which our perturbation results apply. In 
addition, we point out an exact relationship connecting the diagonal part of the Green matrix 
to the asymptotic behavior of solutions. With both of these tools it is possible to identify an 
Agmon metric, in terms of which, in some situations, any decreasing solution must decrease 
exponentially. 

A discussion of the discrete Schrordinger problem and its connection with orthogonal poly- 
nomials on the real line is presented in an Appendix. 

1 Introduction 



In this article we address the asymptotic behavior of solutions to linear difference equations of 
Schrodinger type, as the index n tends to infinity. We prove exponential dichotomy theorems and 
refined approximative expressions for the growing and subdominant (i.e., decaying) solutions, which 
have controlled errors. 

We begin by approaching the subject as a perturbation analysis, showing that if two Schrodinger 
difference equations have potential terms that are sufficiently close, then they are asymptotically 
equivalent in the sense of [12j . that is, there are solution bases for the two problems with the 
same behavior at infinity. The expressions obtained by the perturbation analysis are not merely 
asymptotic, but convergent for large but finite indices n. We follow with a classification of the 
possible asymptotic behaviors and some more estimates, including some cases where the asymptotic 
behavior of solutions does not match that of the comparison equation but can nonetheless be 
characterized. 
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Of course, when faced with one particular equation, comparison theorems are of Hmited use in 
the absence of a good equation to which one can compare. We therefore present some methods for 
constructing such equations after the perturbation analysis. Finally, we present some examples and 
remarks about connections with orthogonal polynomials. 

Let A denote the discrete second-difference operator on the positive integer lattice. We stan- 
dardize the Laplacian such that (A/)„ := fn+i + fn-i — '^fn for / = (/„) S ^^(N), and consider 
pairs of equations of the form 

(_A + F)^ = 0, (1.1) 
(_A + y> = 0, (1.2) 

where the potential-energy functions V, are diagonal operators with real values Vn and 
respectively. (Complex Vn and could be allowed with, for the most part, only straightforward 
complications, but we prefer to keep the exposition focused.) 

Our first aim is to find conditions under which the solutions of (|l.ip have the same asymptotic 
behavior as n — )• oo as those of the comparison (jl.2|) when the potential energies V and Vq are close 
in a suitable sense. One application of the analysis is to the asymptotic behavior of eigenfunctions, 
in which case instead of (jl.ip one could write 

{-A + V - E)tl; = (1.3) 

and (—A + V^ — E)(p = for some real eigenvalue E. Again, for simplicity we shall absorb E into 
the definition of V, with no material restriction, because we consider the full set of solutions to (|l.ip 
without restricting to eigensolutions of a particular realization of —A + y as an operator. Those 
interested in decay properties of eigenfunctions should systematically replace Vn in this article by 

Vn-E. 

We do assume, however, that among the solution set of the comparison equation (jl.2p there is a 
distinguished solution that decreases at infinity, unique up to a multiplicative constant, and we follow 
the nomenclature of ordinary differential equations in referring to such solutions as subdominant. 
(The term recessive is also frequently used.) We recall at this stage that if V has a constant value 
^oo ^ [~4,0], then explicit solutions are easily found, and it emerges that (—A — Voo)4' = has a 
subdominant solution, indeed, one that decreases exponentially (see Example 16. ip . Conversely, if 
V =Voo € [—4,0], then there are no subdominant solutions. The significance of the interval [—4,0] 
is that it is the spectrum of A. 

Remark 1. Equation (II. ip is invariant under the transformation 

^ (-l)>n (1.4) 
Vn ^ -4 - Vn, (1.5) 

as can be easily verified. Because of this, any fact proved under the assumption, for example, that 
Vn > has a counterpart for Vn < —4. We shall use this remark to avoid repetition in some of our 
proofs. 

When (jl.2p has a subdominant solution, it will be denoted (fixing an overall constant), and 
ordinarily we shall identify a second, independent solution as cp'^ . We recall that the Wronskian of 
two solutions of a discrete Schrodinger equation, 

W[<P.,(t>+] ■■= ^n<Pn+l - 4>n+l<t>t, (1-6) 
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is independent of the coordinate n, analogously to a well-known fact for Sturm-Liouville equations. 
(See, e.g., [1].) In terms of difference operators V^, 



fn+i - fn and V /„ : 



fn fn—1- 



(1.7) 



the Wronskian can also be expressed as 



W = cl>-{V^<P+)-<P+{V^cl>-). 



(1.8) 



For future reference we recall some simple relations for the difference operators: 



1. V+(V-/„) = V-(V+/„) = A/„; 



2. (Chain Rule) V+{fg)n = (V+/„)(7n + fn{V+9n) + (V+/„)(V+5n); 



3. (Chain Rule) V~{fg) 



(V-/0<7n + fn{V~gn) - {V^ fn){V- On). 



In comparing (II. ip and ()1.2p the behavior of solutions at infinity will be examined from several 
points of view. First, we study the asymptotic behavior of solutions under the effect of small 
perturbations of the potential as a fixed-point problem. We consider the solutions of (II. 2p as known, 
and use them as the basis for a (convergent) variation-of-constants calculation of the solutions of 
(II. ip . Then we introduce a factorization of the equation satisfied by the coefficients in that scheme 
in order to get a detailed understanding of how they converge. 

Thereafter we present a new and efficient discrete variant of the Liouville-Green (WKB) approx- 
imation [22], so that for a given potential V a comparison equation (|1.2p can be found for which 
the asymptotics are explicitly known, and consequently the behavior at infinity of solutions of (|l.ip 
is explicitly determined, with controlled errors. As an alternative, following [ill ITj 18]. we explore a 
set of related exact relations based on the diagonal of the Green matrix and their consequences for 
the behavior of solutions at infinity. 

Finally, the reader may refer to the Appendix for a discussion of the relation between orthogonal 
polynomials and second-order difference equations. There the connection between ratio asymptotics 
of orthogonal polynomials and the results of Geronimo-Smith [13] will be discussed, and it will 
be shown how solutions of the discrete Schrodinger equation can be represented by orthogonal 
polynomials of the first and second kind. 

We are far from the first to consider these questions, and like other researchers we mimic the 
better-developed theory known for Sturm-Liouville problems. Let us close the Introduction by 
placing our work in the context of the earlier literature. 

A systematic study of certain difference equations dates from Poincare [21]. In his work and in 
that of Birkhoff [1] asymptotic analysis was considered for equations using what would nowadays 
be termed transfer matrices of special types. A rather satisfactory understanding of the effect of 
small perturbations on stability questions for equations using transfer matrices, with dichotomy 
assumptions on their eigenvalues, was developed in [23l [101 E], some of which is recounted in the 
monograph by Agarwal [1], which is a good source for showing how many of the standard facts 
from Sturm-Liouville theory can be ported over to the discrete setting, in particular, the technique 
of variation of constants. Coffman |10j and Benzaid-Lutz ^ studied product solutions, and in that 
regard prefigure in a rough way what we do in Section jH The main results of [3] were discrete 
analogues of Levinson's fundamental lemma [20] for the asymptotic expression of the solution of a 
perturbed linear differential equation. In [3] the authors considered difference equations of the form 




(1.9) 
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Here A{k) is an N x N diagonal matrix with non-zero diagonal entries {^j{k))jL^ that satisfy a 
certain dichotomy condition. They further considered 

x{k + 1) = [Ao + V{k) + R{k)]x{k), (1.10) 

again where Aq is diagonal and V and R satisfy certain bounds. Their results apply widely to 
perturbed difference equations, but not readily to (jl.ip and (jl.2p : As we shall see in (j2.3p . the 
transfer matrices in the present article are of the form 

i+Mn = i+^^^^( (i.n) 



w 

which are neither diagonal nor diagonable ii Vn — 0. In fact, 1 is the only eigenvalue of the 

matrix / + M„ and it has geometric multiplicity one. Moreover, the term {4>n)'^ the lower 

left corner typically diverges as n — ?• oo. 

Trench [281 129] succeeded in giving conditions for the asymptotic equivalence of the solution sets 
of (jl.ip and (|1.2p in the sense considered by Hartman and Wintner [T7], and seems to have been 
the first to realize that a good criterion for equivalence relies on an analysis of the expression 

Jk ■■= cpt^k i^k - Vk') , (1.12) 

(in our notation). In [0], following Trench, a necessary and sufficient condition for asymptotic 
equivalence for some difference equations related to (jl.ip and ()1.2p is spelled out in terms of J. 
Although we bring different methods to bear on asymptotic equivalence in the following sections, 
iP norms of (|1.12p and similar quantities remain central; see Theorems 12.21 13. 2| and 13.31 One 
could interpret these norms as traces of operator perturbations like those occurring in studies of 
spectral-shift functions (e.g., see [H]), leading us to speculate that direct connections between the 
spectral-shift functions and behavior at infinity could be found. 

After a discussion of asymptotic equivalence, we take advantage of the specific Schrodinger form 
of the equation, and construct comparison equations having product solutions of a certain structure, 
inspired by the classical Liouville-Green, or WKB, approximation. Of prior work on discrete versions 
of the Liouville-Green approximation we single out that of Geronimo and Smith [13j . which was 
inspired by some earlier work of Braun [5]. Geronimo and Smith studied a somewhat more general 
equation than (jl.ip . 

dn+lVn+l - qnVn + Vn-l = 0, (1.13) 

where d„ and are sequences of numbers with (i„ / for n = 1, 2, . . . , and pursued a Riccati 
analysis for solutions in product form. In Section U] we identify a more explicit and efficient product 
scheme along the lines of the Liouville-Green approximation as presented in [22], to which we apply 
the perturbation analysis developed in Section [2j Yet another article with Liouville-Green analysis 
using products is [6], in which an explicit semiclassical parameter appears, and the Green matrix is 
studied in a product form and used to prove refined stability results for nonhomogeneous difference 
equations. In the following subsection we relate the discrete Liouville-Green approximation to 
the diagonal of the Green function, following ideas pioneered in [Tl], which have previously been 
somewhat developed in the study of difference equations by Chernyavskaya and Shuster [3 [HI [6] . 

2 Variation of constants and behavior at infinity 

We begin by casting the problem of understanding the asymptotic dependence of solutions at infinity 
as a problem on a certain weighted Banach space, following ideas of [16] in the continuous case. 
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which was in turn inspired by [T7]. Suppose that V is close to another potential such that the 
solutions to (—A + l/'^)(/> = are understood, in the sense that a pair of independent solutions (p^ 
can be identified, including a subdominant solution G A perturbation analysis can be based 
on the following way of connecting the solutions of (jl.ip and (jl.2p . 



Theorem 2.1. Let V and be two potential functions, and let (f>^ be independent solutions to the 
equation ()1.2p . We may represent any if) as a linear combination of (j)^ with variable coefficients 

V'n = On^n +<0n- (2-1) 

Then is a solution to the equation (ll.ip if and only if we may find sequences (a^)5^i that satisfy 
the following two conditions: For all n > 1, 

(V-a+)</.+_i + (V-a-)0;_i = 0, (2.2) 



+ 



w V-«« 

under the convention that 0^ = 0. 

Proof of Theorem \2.1[ Suppose that ip is a solution to (jl.ip . Since the expression ()2.ip has two 
degrees of freedom we have the liberty to impose a second condition on the coefficients to so that 

V-Vn = a+(V-0+) + a~{V-(p-). (2.4) 

Observe that (12. 4p implicitly sets the following expression to zero 

(V-a+)<A+ + iV~a-)cp- - (V-a+)(V~<A+) - iV~a-)iV~<P~) = 0; (2.5) 

by direct expansion (12. 5|) is equivalent to (j2.2p . 

Now compute Ai/;„ = V"''V~'i/'n based on the expression (j2.4p . By the chain rules for V^, 

AV„ = (V+a+)(V-(/>+) + a+(A<A+) + (V+a+)(A</>+) 

+ iV+a-)iV-cP~) + a;(A0;;) + {V+a-){AcP~). (2.6) 

By substituting A^/;„ = Vn.ipn and A(/)^ = V^cp^ into ()2.6p . we obtain 

(Vn - v^)ia+cp+ + a;;<A;^) = (v-<A+ + v;V+)(v+a+) + (v-cp- + y„V;;)(v+a;). (2.7) 



In order for the coefficients of V"'"a„ in (|2.7p and (j2.2p to match, we multiply (j2.7p by (p^ and (j2.2 
by (V~(/)~ + T^c/)"). Then we subtract one from the other and get 

(V+a+)VF = {Vn - K?) (<A+a+ + <Pn = (K - K°)V'n</'n , (2-8 



where W is the Wronskian as defined in (jl.6p . 

Similarly, we match the coefficients of V"^a~ in (j2.7p and (|2.2p by multiplying ()2.7p with (p^ and 
(lO) with (V~(/>+ + (/>+). Then we obtain 



(v+a;;)(-M^) = (Vn - v;0) (<^+</'+a+ + 4>n^t<) = (K - K°)V'n</'+. (2.9) 

Putting (12. 8p and ()2.9p together in matrix form, we arrive at (12. 3p . 
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For the implication in the other direction, suppose that the sequences satisfy ()2.3p and ()2.2p . 
By direct expansion of (j2.ip . we find that 



Atpn = Vni^n + (right side of <^J^) + V+ (left side of ([22])) • (2.10) 

With ([23]), and follow. If we now apply these relations to the right side of ([2T1), it 

becomes 

T/ _ T/O 

[(V-,/.+ + T/V+)Vn0; - {V-K + V^rn)i^nCt>t] , (2.11) 

which simplifies to 

i^n^^^^[(t^-V-<p-^ - 'P-^'^-^n] = (K - V^)i^n. (2.12) 

Returning to (|2.10p . we conclude that A-i/'n = Vni^n for all n only if V"'' {{"^^ CL^)(t>n~i + ^n )'Pn-i) 
for all n, or in other words, when the expression on the left side of (j2.2p is a constant. Finally, 
note that since its value is zero when n = 1, it is zero for all n. □ 

2.1 Convergence of in a suitable Banach space 

An advantage of the variation-of-constants approach to asymptotic equivalence over the methods 
of |28l [28} [9] is that it provides a rapidly convergent iterative scheme with error estimates that can 
be made explicit. To set it up, we introduce the notation 

at\ , . Vn-V^ 



^_^and /3. = ^V^, (2.13) 
where W is the Wronskian as in (jl.6p . We shall regard a as an element of a weighted Banach space, 



X 



BN:={y^ = (£.]■■ ||X|U := sup (|(<A+)'X+| + |X-|) < oo ^ . (2.14) 



n>N 



Substituting the expression for il) into (|2.3p . we calculate 



or, by summing (j2.15p . 



(V-a„) = (2.15) 



a.-a.+i ^M_(^+)2 _^+^-)^" 



— 



Formally letting I ^ oo, an+e — )• , we therefore define a linear operator A4 by 
The convergence of the coefficients in the Banach space proceeds as follows: 
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Theorem 2.2. Suppose that ()1.2p has a solution basis (p^ such that lim.n^ao4'n ~ 0' I't'nl ^-^ 
monotonically nondecreasing for sufficiently large n, and f3n (cf. (|2.13p ) satisfies + IfAn*^?! P) ^ 
Then for N sufficiently large, A4 is a contraction on Bj\f. Consequently, there exists a unique 
solution of (jl.ip such that 

where lim„_s.oo «n = and lim„_s.oo = 1- Moreover, if we define 'tp' := maxm>n |</'n I; then 

i^n = 4>n + rn^n , (2-17) 

with lim„^oo rn = 0- 

Remark 2. If |(/)~| is monotone nonincr easing, then we may simply write ip~ = (1 + rn)^'^^ with 
hm„^oo fn = 0. If the product (pn4>n bounded, it suffices for this theorem to assume that /3„ G 
circumstances under which this is guaranteed are discussed below in Lemma [2.211 Theorem 14.31 and 
Theorem 14.21 

Proof. For N sufficiently large, we claim that M is a strict contraction on B^. To see this, 

for sup„>;v(|X+| + |X„ I), and observe that ||X||gjy = 



we introduce the shorthand sup 



X„ 



SUPm>Af 



l'^" 



Thus 



||-A/^X||BJ^ = sup 



n>N 



sup 

n>N 



1 j 2^ ^"+^' 

^ k=0 

oo 
k=0 



K+k'l^n+k 



X 



n+fc 



+ P 
1 



n+k^^n+k 



1 



-1 



sup 

n>Ar 



fc=0 



n+fc 



^n+k^n+k 



< W^Wbm sup 

n>7V 



l' 



fc=0 



<2||X||B^^|/3;v+fe|(l + 



''N+k'f'N+k\ ) 



fe=0 



We then ask whether there is a solution to 



3-n, 



(Ala) 



and conclude by the contraction mapping theorem that there is, for N large enough that 

oo ^ 
^WN+k\{l + 



'^N+k't'N+kl'^) < n- 



k=0 
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Indeed, therefore (1 + A^)a = is uniquely solved by the norm-convergent Neumann series 

-=(f:i-^A(^)- (2.18) 



Being dominated by a geometric series, the convergence of (j2.18p is exponentially fast. 
For the final statement, we need a lemma about the operator M: 

Lemma 2.3. Suppose that \x^\ < Ci|(/)^p and \x~\ < C2. Then 



and 



where M := ||/3„(1 + |(/.+</)-| 



\Mx\+ <M {Ci + C2M~f (2.19) 
\Mx\- <M {C1 + C2), (2.20) 



The lemma is an easy estimate from the definition of M. 

The proof of the final statement of the theorem then requires choosing N sufficiently large that 
the coefficients in the conclusions of the lemma are small enough that M is a contraction, and 
summing the Neumann series ()2.18p . □ 

2.2 Construction of a second solution and estimates of the product of the two 
solutions 

Since understanding the asymptotic behavior of the solutions of the perturbed equation requires 
knowledge of a full set of independent solutions </'n} to the original equation, we recall a 
standard reduction-of-order formula showing that the subdominant solution determines a second, 
independent solution, which grows at infinity. (E.g., the text |17j treats this argument in the 
continuous case in §XI.2, and it can be found in the discrete literature in numerous places, including 
|28t [9].) The following simple formula does not require a subdominant solution, only one that is 
nonvanishing. It is true by direct verification that ip^ solves (jl.ip and that W = 1 when n = m. 
(Of course it is derived by positing that ip^ = ^ substituting, and using the Wronski identity 
to determine 7,1.) 

Lemma 2.4. (Standard) Suppose that (jl.ip has a solution that is nonzero for all n, n = m, . . . , M . 
If 

if n = m 

V'+:= 1 ^ (2.21) 

't'nZ^-^-^ ^fn>m. 

then tp^ is an independent solution of p.ip on the interval [m,M]. 

The lemma has some simple but useful consequences: 

Corollary 2.5. Suppose that for some a, \a\ > 1, (jl.ip has a solution such that a^ip~ 1 as 
n — )• 00. Then 

• Every solution ipn of p.ip that is independent ofip~ is exponentially increasing, i.e., a~'"'ipn — ^ 
C 7^ as n —)■ 00. 
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• For any solution ipn, the product ipni^n bounded independently of n. 

• Given any boundary condition of the form atpi + 6'i/'2 = 0, ifO^ sp(— A + V), then the Green 
matrix for (—A + V) on n > 1 is uniformly bounded. 

Proof. Because every solution to (jl.ip is a linear combination of and V'n defined in p.2ip , it 
suffices to show the first two statements for ipn = V'n ' which behaves asymptotically like 

n-l 



•^a^^+i. (2.22) 



k 

Since this geometric series can be bounded above and below by (Ci + C2 a"" J" ^ a'^^'^^dx = 
^a"-i + 0(1)) Ci + C20" the first two statements follow. 
If 7^ sp(— A + V), then the Green matrix is defined, and 

^ ^min(m,n)^max(m,n) /„ „„s 

where satisfies the boundary condition at n = 1,2. Since this is bounded on any finite set of 
indices m, n, the third statement follows from the asymptotic estimate of "0™ the first statement. 

□ 

An important case where these estimates apply is captured in the following. 

Corollary 2.6. Suppose that for some constant Voo ^ [—4,0], V — Voo G ■ Then there is a solution 
to (II. ip of the type ~ a~" for \a\ > 1 and an independent solution ~ a", and the statements 
of Corollary 12.51 apply. 

Proof. We can apply Theorem l2.2l and Corollary 12. 5l once it is observed that the comparison equation 

(-A + yoo)0n =0 



has an exponentially decreasing solution, viz., assuming Voo > 0, = a " for a = ^ yVoo + \/V^ + 4 Vex 
(The case Voo < —4 similarly has an exponentially decreasing solution, according to Remark [TJ) 

□ 

Further conditions for the existence of exponentially decreasing and exponentially increasing 
solutions may be found in |32j . 

The existence of a more rapidly decreasing solution has similar implications: 

Corollary 2.7. Suppose that for some a > 1,6 > 1, equation (|l.ip has a solution such that 
a^^i^n ~^ ^ as n ^ 00. Then 

• Every solution tpn of (jl.ip that is independent of tp' increases rapidly as n — t- 00, but 
o.~^^il^n — )• as n — )■ 00. 

• For any solution ipn, the product ipnipn bounded independently of n. 

• Given any boundary condition of the form atpi + bip2 = 0, i/0 / sp(— A + V), then the Green 
matrix for {—A + V) on n > 1 is uniformly bounded. 
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The proof of this is similar to that of Corollary 12.51 but details will be left to the interested 
reader. On the other hand, subdominant solutions that decrease only polynomially fast do not lead 
to as strong control of the products or of the Green matrix: 

Corollary 2.8. Suppose that for some a > Equation (jl.ip has a solution such that n^ipn ^ cls 
n — )■ oo. Then for any solution ip that is independent of , \'tp~tpn\ ~ Cn. 



Proof sketch. The argument being familiar from the proof of Corollary 12. 5^ we content ourselves 
with the application of Formula (|2.2ip . As before, we may as well assume that ■0n = V'n defined 
by that formula, which is of the form 



n-l 

— n 

n 

k 



""^k^ik + lf. (2.24) 



This is asymptotic to n~°' ^ x'^"'dx = 2a+i ~^ o{n°'~^^). The claimed estimate for the product 
results when this is multiplied by ip~ . □ 
Some converse implications, by which the boundedness of 4>n4'n controls the asymptotic behavior 
of solutions, will appear in Section HI 

3 Refined asymptotic estimates 

Let ip be a non-trivial solution to the equation —A + V = 0, when — is small, and ^pn = 
^n4'n ~^^n4'n- ^^is section, we provide a classification of the parameters and a~ and describe 
how they converge. 

To begin, we prove a preliminary classification of a"*" and a~ , distinguishing the exceptional cases 
where the perturbed solutions only depend on one of the comparison solutions in ()2.ip for large n: 

Proposition 3.1 (primary classification). Suppose I3n<pn4'n ~^ 0- Then for any non-trivial solution 
ijjn = a+(/)+ + a~(t)~ with a^ satisfying (j2.3p . one of the following must be true: 

1. Given any integer N , there is an integer p > N such that both and a~ are non-zero (this 
will be treated in Theorem \3.2\ below). 



2. There is an integer pQ such that 



aJ+, = a+^0, a-+, = V/c > 0. (3.1) 

3. There is an integer pQ such that 

a-+, = a-^0, a++fe = VA; > 0. (3.2) 

Proof of Proposition [3Tl Recall the definition in ()2.3p . Observe that detMp = TrMp = 0. Hence, 
for all p G N, 

det(/ + Mp) = 1 + TrMp + det Mp = 1, (3.3) 

which implies that 

/^+\ 

= for some m <^ a^ = 0, Vp G N <^ ai = 0. (3.4) 
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Since ip is not the trivial solution, for any p either 7^ or 7^ (or both). 

Suppose we are not in Case 1. Then there exists Nq such that for all k > Nq, either = or 
= 0. Without loss of generality, suppose for some large p > A^'o, a+ / and 1 + /3;j(/)^0^ / for 

all k > p (this is possible because l^k'Pk^k ~^ *-*)• ^^^^ show that this corresponds to Case 2. 
Observe that 

which implies Op+i 7^ and as a result, ap_^_-^ = 0. Apply the same argument recursively to obtain 
(|XTD □ 

Assumptions 3.1. Now we focus on Case 1 of Proposition 13. 1[ Without loss of generality, we 
assume af,a]^ / and sup„ |/3n^+(/>~| < 1. The latter assumption, together with the assumption in 
Theorem 13.21 that Yl'^=i \l^n4't'Pn \ < guarantees that 7^ for all n and that lim„_s.oo / 0. 



A key observation here is that the recurrence matrix I + M^ in (12. 3p can not always be diagonal- 
ized, making it impossible to utilize existing techniques in the perturbation theory literature, which 
heavily relies on the fact that the transfer matrix can be diagonalized (see, e.g., Benzaid-Lutz p] 
and the discussion in the Introduction). 

A key observation is that the recurrence matrix can be decomposed into the sum of a lower 
triangular matrix and an upper triangular error matrix: 

=: Gn + En. (3.7) 
The advantage of such a decomposition is as follows: let S„ be defined recursively as 



(3.8) 



n-l 

-M<Pif 11(1 + Pj<P^(t>j) + (1 - Mi<p-)^n-i, n>2- 

_ -/3i((/>+)2, n = l, 

and let 

n n 

n± = n(i±/^^'^/'^7) ■= Upf (3.9) 

i=i i=i 
Under such definitions, we have a closed form for the product GnGn~i ■ ■ - Gi: 

GnGn-l • • • Gi = G„ ( "^^ ) = ( rr- ) • (3.10) 

While Gn is an approximation for the recurrence matrix / + we shall show that the product 

GnGn~i • • • Gi will serve as an approximation to the actual recurrence relation (|2.3|) for under 
Trench- type conditions on Vn — Vn- 

Theorem 3.2. Let (jx^ be independent solutions to the difference equation —/S. + V^ = 0. Consider 
a potential V such that 

oo 

sup < oo and \J |(14 — ^)(<An'An)| < oo. (3.11) 

n 

n 

If either ^1^=1 I'/'n P < or ^1^=1 ~ '^n\ < oo, then any non-trivial solution ip to the 
equation —A -1-^ = can be written as tpn = o-n^'n + ^n'Pn such that one of the following is true: 
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1. There is an integer p such that 

«P+A: = </0, a;+, = Vfc>0. (3.12) 

2. There is an integer p such that 

a;+, = ap/0, a+^, = Vfe > 0. (3.13) 

3. There exists a constant foo ^ such that = H-nfoo + o{l) and a~ = Sn/oo + o(l). 

4- There exist constants f^,f^ with / such that a+ = n+/+ + o(l) and a~ = ^nfto + 

/oo + 0(l)- 

In Example 16. 11 we will construct a potential V such that sup < oo but S„ fluctuates as n 

n 

goes to infinity. 

Next, we prove a result such that may go to infinity: 

Theorem 3.3. Let the definitions be the same as in Theorem \3.2\ and V be a potential such that 

oo 

Vn - K? > {or < 0) Vn G N and ^ [(K - T/°)(0+O| < oo. (3.14) 

n 

Moreover, we require that sup^ |(</>~)^S„| < oo. 

If lim = oo, then exactly one of the following must be true: as n ^ oo, 

n— >oo 

1. a+ = n+(a+ + o(l)) / and a;^ = Sn-i(ai + o(l)); 
£ — )■ and a~ — )■ aj^ ^ 0. 
//sup|S„| < oo, the reader may refer to Theorem \3.2[ 



Proof of Theorem \3.3[ Let f^+i and f^+i be defined implicitly in (j3.15p below: 



n+1 



(3.15) 



First, we want to prove that 

sup|/+| + |/-| <oo. (3.16) 

n 

Then we prove that exactly one of the following must be true: 

1. lim /+=:/+ / and lim /" = 0. 

n— >-oo n— >oo 

2. lim fn=-f^ exists and lim /~ / 0. 

n— >-oo n— >oo 

To begin, we observe that 
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which, by the triangle inequahty, imphes that 



/3„(0-)2n;_i 



n+ 



(3.18) 



Jn+1 - in — TT^ TT^ — — • (O.iyj 



n- 



n-l 



By (1131), 
Thus, (|3.19p becomes 

/n+l ~ /n 

which imphes that 



a„+i - Pn a„ = -/3n(0+)^a;t = -/3„((/)^ )^n+_i/+ = (S„ - p„ S„„i)/+. (3.20) 



s„(/+ - /++i) _ ^.(<^-)2s„(E„_i/+ + n-_j- 



|/„"+il<|/„1 1 + 



/3n(0;^)2s, 



+ 1/^ 



(3.21) 



(3.22) 



Add p.lSp to (j3.22p . Since Y^=\ \(^n'Pt4'n I < °o, 11^ converges to a non-zero hmit and — >■ 1. 
Moreover, sup^^ < oo, so there is a constant K such that 

I/++1I + l/'+il < (1 + K\(3n{cp-f\) + |/-|) , (3.23) 

which imphes (|3.16p . 

Dichotomy 3.1. There are only two mutually exclusive possibilities for and /~: 

1. For any pair consisting of an integer N and a constant M > 0, there exists an integer p = 
p{N, M)> N such that 

M\f;\>\f-\. (3.24) 



2. There exist an integer A'^o and a constant Mq such that 

|/+|<Mo|/-| V?i>7Vo. 



(3.25) 



Suppose we are in Case 1. Note that (|3.24p implies that /+ 7^ 0, because if = 0, then 
fp=0 which implies = a~ = 0, and hence = 0. This contradicts the assumption that ip is 
a non-trivial solution. 

Let p be the integer given in (|3.24p . We shall specify the choice of and M later in the proof. 
Let 

vn = (3.26) 

Jn 



Note that by the triangle inequality and (|3.17p . 



fp+i 



> 1 



f+ _ f+ 



>l-Ki|/3p(</>;)2|(l + |rp|) >0. 



(3.27) 
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Thus, 7^ implies /^^^ 7^ 0. Furthermore, by inverting (j3.27p . we obtain 



p+1 



< 



l-Ki%{4>pY{l + \rp\ 



< Ko. 



(3.28) 



Clearly, both and rp+i are well defined as fp,fp^i 7^ 0. Observe that by the triangle 
inequality, 



kp+i - rp\ < 

By by dMID, 

fp+i ~ fp 



fp+i fp 



J V 



p+1 



+ 



/p+1 fp 



f? 



fp 



f^ 



Jp 



fr 



p+1 



(3.29) 



fp^ 



< 



Similarly, by ()3.17p . 

/p+1 ~ fp 



fp^ 



< 



n+ 



n 



Sp_i| + |rp||n;_J < (l + lrpDKal/SpCcA;)'!- (3.30) 



+ \rp\\U^_^\) < K4(3p{cp-f\{l + \rp\). (3.31) 



By the triangle inequality, 

kp+i| < kp+i - rp\ + \rp\ < \rp\ 



1 + 



K^il + \rp\)\f]p{ 



l-Ki\l3p{^p)^{l + \rp 



+ 



l-Ki|/3p(0p)2(l + |rj, 



^1 (l + M + |r,p)ir5lWp)^l 

'"""'^ l-Ki|/3p(0p-)2(l + |r,|)| ■ ^ ■ ' 

In particular, if |rp| < 1, then 1 + \rp\ + |rpp < 1 + 2|rp|. Besides, when p is large (which will be 
the case). 



1 - i^i|/3p(</>-)2(l + |rp|)| > 1 - 2Ki\Ppicl)- f\ > 1/2. 



Hence, p.32p becomes 

|rp+i| < \rp\ + (1 + 2\rp\)Ke\(3p{4>;f\ = \rp\ (l + 2Ke\Pp{cP;f\) + KeWp{(t>;)^\ 
Hence, by an inductive argument we can prove that 



'k-i 

rp+k\ < I n + 



j=0 



where 



Since either /3fc or 10;, P are summable, ?7fc G Hence, 

00 00 

Pp - n (1 + ^p+j) < n (1 + := < 00 

j=0 j=0 



(3.33) 
(3.34) 

(3.35) 
(3.36) 

(3.37) 
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and 



k-l 



Sp := ^ r]p+j < Kj'^ r]p+j as p ^ oo. 

Here is how we choose M and N: given e > 0, choose M, N such that 

\rp\ < M < and sup Sk < 

^i^oo k>N ^-foo 

It is guaranteed that there exists an integer p > N such that \rp\ < M. By p.35p . 

kp+fcl < -Poo (kpl + Sp) < Po> 

In other words, if (j3.24p is true, then lim r„ = 0. Apply this to (13.3ip . we get 

-f n 



(3.38) 



(3.39) 



(3.40) 



'n+l 



fit 



<Kr\(3n{<l>-f\(^i'. 



(3.41) 



Since logz is analytic near z = 1, in a neighborhood of 1 there is a constant Kg, arbitrarily close to 
1, such that 

|logz| = |logz - logl| < Ks\z - 1|. (3.42) 

Put z = By 



log 



fj 



n+l 



(3.43) 



Moreover, by the argument following p.27p . we know that there exists an integer p such that fn¥^0 
for all n > p. Therefore, 



J=P 



f 



That implies the existence of lim„_j.oo log/+, and 

lim /+ := /+ / 0. 

Together with the proven fact that r„ — t- 0, we obtain 

lim fn = 0- 



(3.44) 



(3.45) 



(3.46) 



Next, suppose (j3.25p is true. If = 0, then = for all n > Nq, which contradicts with the 
assumption that -0 is a non-trivial solution. 

Now suppose there is an m > Nq such that /~ ^ 0. Divide both sides of (j3.2ip by /~. Then 
we obtain: 

fm+l 



fn 



1 



(3.47) 



Using the same argument as in Case 1, we can prove that / implies f„^_^i / 0, which 
allows us to apply the same logarithmic argument to prove that 



lim /„ := / 0. 



(3.48) 
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By (I337D, 



= ^1 + /3„(<^-)^^) /+ + yM<l>nf^j fn = (1 + 0(|/3n(<A;)'l))/+ + 0( l/^n (0 | ) • 

(3.49) 

Hence, 

lim fn = /(i exists. (3.50) 

□ 

Proof of Theorem 13. 3[ The proof is very similar to the one of Theorem 13. 2i The only difference in 
the proof is that instead of f^ and r„ we consider 

fn ■■= 4^ and r„ = (3.51) 
In place of |r„+i — t„| in (13.29p . we replace it with 

(3.52) 



and make use of the fact that < 1 for aU n. □ 

4 Construction of Comparison Equations 

In this section we turn to the problem of determining the asymptotic behavior of solutions of (jl.ip 
as n — 7- oo given a potential Vn, where Vn can be either bounded or unbounded. We shall construct 
explicit comparison equations with respect to which we can call upon the perturbation results of 
the earlier sections of this article. The construction will require a discrete replacement for the 
Liouville- Green (familiarly, WKB) approximation, which is a well-known and quite useful tool for 
this purpose in the setting of ordinary differential equations |17^ 122]. 

Our ansatz is that given an equation of the type (II. ip . a related equation is to be sought for 
which the solutions are of the form 

n 

<t>t=^n\{Sf\ (4.1) 

1=1 

Recall that in the Liouville-Green approximation to ordinary differential equations of Schrodinger 
type a comparison is made to a similar equation having a solution basis in the form 



|22j . In common with previous authors, we replace the exponential function containing an "action 
integral" by a product of the quantities we designate Sn, but we innovate with an additional prefactor 
z„, to be specified below in (14. Sp . This is designed to bring simplifications in the discrete case 
analogous to those resulting from the prefactor V{x)~^/^ in the continuous case. 

Before we state the main results of this section, we pause to point out a connection between the 
factor Zn and the Green matrix for the Schrodinger operator —A. + V, viz., for n > m, 

n ^ 

Gnm = <Amin(m,n)'^max(m,n) ^ 11 ' ('^•^) 

t=m+l 
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for which (—A + V)G is the identity operator, by a direct computation. In case n = m, 

Gnn = ^n- (4-3) 

In the following section we study the diagonal elements of the Green matrix and show that they 
are directly related to the behavior at infinity of solutions and to the notion of an Agmon metric, 
(cf. [m Section 4]). 

In order to determine z^, we recall the constancy of the Wronskian of solutions to equations of 
the type (jl.ip . To simplify the discussion, we take W = 1, which can be arranged by scaling. Given 
our assumptions it implies that 

ZnZn+l (Sn+1 - -J—] = 1- (4.4) 



Sn+1 

Guided by the case of a constant potential, we expect that if Vn is well-behaved, then a good 
choice for Sn is one of the solutions of Sn + = Vn + 2. This turns out to be adequate in some 
bounded cases, but a more sophisticated choice is necessary when, for example, Vn is allowed to be 
unbounded. We remark that the choice is not unique, because different choices lead to the same 
asymptotic behavior if the comparison potentials they lead to are sufficiently close. Our discussion 
will proceed under the supposition that > for large n; the case where Vn < —4 for large n is 
similar with the systematic sign changes mentioned in Remark [TJ 

To determine the best choices for Sn and z„, we consider the equation 

Sn+^= bn, (4.5) 

which is effectively a quadratic, and let 5„ be the root of larger magnitude, i.e., 



2 

We observe that the relationship 



Sn = ^ , where |5„| > 2. (4.6) 



1 



Sn-^ = (4.7) 



necessarily follows. 

To be consistent with the Wronski identity (j4.4p we must set 



4)(6L,-4)... 



for all n, where the constant Cz will be chosen below. (We clarify that the prefactor simply alternates 
between Cz and its reciprocal, depending on whether n is even or odd.) 

The comparison functions (f)^ both solve a Schrodinger equation with potential Vn given by 

K:=^ = ^5.^i + ^^-2. (4.9) 

(This equation is true by direct substitution for (f)~^; to see that is it also true for (f>~ requires also 
substituting from (j4.4p : cf. a similar argument for (jS.Sp in ^) We shall in fact show that there is 
a choice of ways to choose bn that will lead to a sufficient convergence rate oi Vn — Vn, and that 



17 



the the logarithm of the quantity S'„ can be regarded as an Agmon metric [21 [18] controlling the 
behavior of solutions (p of (jl.ip at infinity. 

For clarity, we first consider the case where Vn is bounded and > C > for all n > Nq. 
Without loss of generality we may assume that A'o = 1, because this does not affect the large-n 
behavior of a solution basis. (This simply allows us to avoid choosing phases for some square-roots 
of quantities that might otherwise not be positive.) 

In the case of bounded, slowly varying potentials Vn, Theorem 14.11 contains estimates for 5„ and 
Zn and uses them to control the solutions and Green matrix of the comparison equation {—A+V)(j) = 
0. The construction in Theorem 14.11 is guided by the special case of a constant potential. 

In Theorem 14 . 2 1 that follows, we shall present a more general result which covers potentials that 
are convergent to a finite limit under more relaxed assumptions on Vn- 

Finally, in Theorem 14.31 we show that the method proposed in this section also works for un- 
bounded potentials that possibly fluctuate. 

Theorem 4.1. (bounded and slowly varying potential) Suppose that for some C > 0, C < Vn £ i°° , 
and that n{Vn+i — Vn) G ■ Choose 

bn = b'n^'' ■.= Vn + 2. (4.10) 

This implies (with a short calculation) that 

Sn = Si^'' := i + 2 + v'K(K + 4)) . (4.11) 

The factor Zn is determined by (|4.8p . Then 
(a) 

S^n^i - S^^^'^ e f. (4.12) 

(b) Vn converges to a nonzero limit Voo as n ^ oo, and 



m=i V ^'^rn-l{V2m-l + 4j 

is well defined through a finite convergent product. 

(c) Under this definition of Cz and the one of Zn in ()4.8p . ZnZn+i = l/\/Ki(Ki + 4) and 

1 



[^oo(v;o + 4)] 



1/4 



E l^. (4.14) 



(d) The comparison potential Vn defined in (14. Oh . satisfies lim„_5.oo Vn = Voo, o-nd the Green matrix 
for —/S. + V and the product (pn't'n uniformly bounded. 

(^) Vn — Vn G and therefore, identifying V with the comparison potential V^ in (jl.2p . the solutions 
of (II. ip and ()1.2p are asymptotically equivalent in the sense of Theorem \2.S\ . 



Proof. The proof for (a) is a direct application of Taylor's Theorem: Following (|4.6p . we consider 
the function f{x) = l/2(x + 2 + ^ x{x -|- 4)), which is differentiable for all x > 0. In particular, if 
x,y > C > 0, /(y) = fix) + R{x,y), where R{x,y) = [y — x){r + 2) r{r + 4) for some r between 
x,y, implying that R{x,y) is uniformly bounded in x,y if x,y > C > 0. Since S^^i = f{Vn+i) we 
can write 

= 5^ + R{Vn, K+l) = S'n'''' + 0(|K+1 - Vn\), (4.15) 
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which proves (a). 

The fact that if the differences Vn+i — Vn are summable, then Vn has a hmit is immediate. To 
estabhsh the convergence of the product (|4.13p . let 5m ■= Vm — Kn-i- Then 



Vm{Vm + 4) = Vm-i{Vm^i + 4) + 6m{2Vm~i + 6m. + 4) 
and since < C < < + 4 for ah n, 

2Vm~i + (5™, + 4 



As a result, 



"Vm-l{Vm-l+^) 
Vm{Vm + 4) 



0{\5m\). 



Vm-l{Vm-l+^) 



1 + 



Vm - Vr 



m—1 



m—1 



(4.16) 



(4.17) 



(4.18) 



Vm(Vm+4) 



^rj^ S (-^ ■ By taking the logarithm in (j4.13|) . the product therefore 



which imphes that In -rj 

Vm-l(Vm 

converges, and is easily seen to be nonzero. This proves (b). 

The same argument for (b) establishes the convergence as m — >■ oo of Z2m and of Z2m+i-, sepa- 
rately. The choice of the prefactor in (j4.13p ensures that the two limits are the same. The more 
precise statement (I4.14p is where the assumption that not only Vn+i — Vn but also n{Vn+i — Vn) G 
is needed. From the definition of Zn it can be seen (by taking logs and using Taylor's theorem) 
that \zn+2 — Zn\ is dominated by a constant times \Vn+i — Vn\ + \ Vn — Vn~i\- Thus — -Zool 

is dominated by a constant times 

oo oo 
m=l k=m 

which by reversing the order of summation equals 



Y,{k-l)\Vk-Vk-i\ < 



oo 



k=l 



by assumption. The other statements in (c) follow by (14. 2p and ()4.3p . 
Finally, by (gj]), 



Rn '■— Vn — Vn — Vn - 

With the aid of gH]), 



Zn+l>->n+l _|_ Zn-l 



y qbdd 



{Vn + 2) 



obt 
-l^n 



bdd 

+1 _|_ Zn-1 



y qbdd 



(4.19) 



Rn — {Sn+l — Sn) + 



Zn+1 



1 Sn+1 + 



Zn-1 



Sn 



Since Sn and Zn both have finite nonzero limits and 5„ — limfc_>.oo Sk and — Zoo are both each 
of these three terms is easily seen to belong to 

□ 

Remark 3. The quantity Cz is analogous to the exponential of an action integral in the continuous 
situation, which shows up in "tunneling" effects. We summarize that in the case where < C < 
Vn and Vn+i — Vn € there is a Liouville-Green basis of comparison functions for (jl.ip . and 
The perturbation method of ^lets the solutions {V'n} determined from that basis through an 
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iteration that converges for all n > for some finite A^. To collect the details in one formula, the 
Liouville-Green basis is of the explicit form 



+ 4)K-3(K-3 + 4) . . . " / 14 + 2 + y^VkiVk + 4) 



±1 



Kx(K + 4)K-2(K-2+4), 



n 



(4.20) 



(dropping the normalization factors or, resp., l/C^). 

The next result concerns potentials that are convergent to a finite limit under more relaxed 
assumptions on Vn- 

Theorem 4.2. (general bounded potential) Let Sn+l/Sn = bn- where bn is a bounded function of the 
potential V such that for all n, Ci > bn > C2 > 2 for some constants Ci,C2- If^n\^n+i — bn\ < 00, 
then 



(a) 
(b) 



+1 



0{\h 



'n+1 



P. 



M 



M 

m=l ^ ^ 



bn\). 



(4.21) 



is well defined through a finite convergent product. If bn 7^ fen+i for all n, then Poo 0. 

(c) If Vn+i — Vn ^ , we can choose bn in a way such that bn+i — bn S and that Vn — Vn ^ . 
Some appropriate choice will be shown explicitly in the proof below. 



Proof of Theorem \4-^ The proof for (a) is a direct application of Taylor's Theorem. Following 
(14. 6p . we consider the function f{x) = 1/2 (x + x"^ — 4), which is differentiable for all x^ > 4. In 



particular, if x, y > 2, f{y) = f{x) + {y — x)R{x, y), where R{x, y) = r j^yjr'^ — 4 for some r between 
X,?/, implying that R{x,y) is uniformly bounded in x,y i( x,y > 2. Since Sn+i = /(6n+i) we can 
write 

Sn+l = Sn + {bn+1 " bn)R{bn, 6n+l) = 5„ + 0(|6„+i - 6„|), (4.22) 

which proves (a). 

Note that (j4.30p can be expressed as 



Vn-Vn = ^ 



bn+1 



Jbl-,1 - 4 



{bl-4){bl_,-4).- 
(6Li-4)(62_3_4). 



(I) 



bn-i 



bl+1 - 4 



(fc^-4)(6L2-4)--- 
(6Li-4)(62_3-4).. 



(K + 2). (4.23) 



{II) 

Clearly, (I) is summable in n because both Sn and bn are of bounded variation. Hence, we are 
left with 

(in - ibl-4)(bl_,-4)... . .4 24) 
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To see what possible choices of 5„ that will give us the desired convergence, we let 

(6Li-4)(62_3-4) 



= . 2" _ \. • (4.25) 



Then we have 
which implies that 



6^+1-4 = J„+iJ„ (4.26) 



'Jn+lJn + 4 _|_ 2^ _ /1 1 Jn — Jri' 

Jn+lJn y 

Therefore, a natural choice for J„ is 



(II) = {Vn + 2) = Jl+ ^ -^"'^ V(-/n+l - Jn)Jn + + 4 - (K + 2). (4.27) 

V Jn+lJn V -'n+l 



+ 4 = (K + 2)' (4.28) 

(or equivalents, J„ = VK(K + 4)), then J„+i - J„ = 0(|K+i - K|) and by (I4.26p . — bn is 
also Od^n+i — Ki|) • Under this particular choice of Jn, 



Wl + , "^'^ ' ^(Jn+1 - Jn) Jn + + 4 = (1 + 0(K+1 - K))(K + 2). (4.29) 

In fact, there are a number of choices of Jn that we can choose from. By (|4.23p above, bn+i — hn € 
and Jn+i — Jn ^ f-^ are are sufficient conditions for Vn — Vn ^ i^- 
We provide a few examples here for the interested reader: 

(a) (geometric mean) Let + 4 = {Vn+i + 2)(y„ + 2). Clearly, under this choice, + 4 = (T^ + 

2)^ + 0{\Vn+l-Vn\) and by (g^S]), 6n+l-6n = 0{\Jn+l-Jn-l\) = 0(| K+1 " K | + | K " K-1 1 ) G • 

(b) (arithmetic mean) Let + 4 = [(K + 2)^ + {Vn-i + 2)2]/2. 

(c) (skipping some K's) For /c > 0, let + 4 = J^,^^ + 4 = (y2fc + 2)^. Then for all n, 
J„+i - J„ = 0(|K+i - K-il) G Hence, bk+i - hk = 0{Vk+i - S 

□ 

We now turn to the case where Vn is not bounded. Here we find that the Liouville- Green 
approximation for the unbounded case simply requires replacing Vn + 2 by the geometric mean of 
Vn + 2 and Vn-i + 2. That is, the canonical choice in the unbounded case is 

1 



Sn-^ = x/(K + 2)(K_i + 2), 
»-'■« 

which, we remark, is equivalent to + 2 when Vn is bounded and slowly varying. The argument 
establishing the accuracy of the Liouville-Green approximation runs much as in the simpler, more 
restricted case, but with correspondingly more complicated details. As before, this choice of Sn is 
convenient but not unique. 

In terms of a general 6„ and the Sn related to it according to ()4.5p . the comparison potential 
(14. 9p becomes 



Vn-Vn 



62^1-4 Sn^/bl - 4 



4)(^n-2 4)--- _(y^^2). (4.30) 



(&Li-4)(62_3-4) 
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Note that by (jilB]) and the fact that 1 / Sn = 6„ - 5„ 



'n+l 



5n+l 



K+1 - 4 



2^/^^+1-4 



1 

+ 2 



SnVW^ 2^6^ 2' 



(4.31) 
(4.32) 



Hence, (j4.30p becomes 



Vn-Vn 



'Jn+1 



+ 



:6Li-4)(62_3-4)... 



iVn + 2). (4.33) 



It turns out that the convenient choice of Cz in a situation where Vn is unbounded is simply 
= 1. The theorem reads as follows: 

Theorem 4.3. (unbounded potential) Let V be a potential that satisfies 



Yl ..1/2 ( ..3/2 + ..3/2 ) < 

Zei = 1 and bn = Sn + ^/Sn > 6e chosen such that 

1 



•^n - ^ = V&n-^ = v'(K + 2)(K-l+2). 



Then 



(a) bn = v'(K + 2)(K-i + 2) + 4, z„ 



1 



VK+2 



and 



^ _ V(K + 2)(K-i + 2) + + (K + 2)(K-i + 2) ^ ^ 



(4.34) 



(4.35) 



(4.36) 



(b) Vn-Vn = 



1 1 

+ 



..3/2 ..1/2 ..1/2.. 3/2 



(c) The Green matrix Gm,n is uniformly bounded. 



Remark 4. The condition (j4.34p is satisfied by unbounded potentials, including some that fluctuate. 
An example is given ^ 



Proof of Theorem (j4.3p . Given this choice of 6^ — 4, we have 

{bl-4){bl_,-A)--- 



(?^Li-4)(62_3-4). 

By the definition of z„ in ()4.8p . this implies (a). 
Therefore, by (|4.33p above, 



Vn + 2. 



(4.37) 



Vn-Vn 



Vn + 2 



'Jn+1 



+ 



bl+1 - 4 



(4.38) 



22 



We apply the relation ^/a — \fh = {a — h)/ {y/a + \/h) to 



bl+i - 4 - 4 ^62^1 - 4 + ^62^1 - 4j - 4 + ^^^^ - 4 

= , , (4.39) 

+ 2)(K + 2) [{Vn+i + 2)(K + 2) + 4] + + 2)(K + 2) 

which is in the order of 0((K+i + 2)^/2 (y^ + 2 )^/^). 

Next, we obtain a similar formula for h^j — 4 and show that it is in the order of 0((V^ + 
2)3/2(T4_i + 2)3/2). Canceling the K + 2 term in (OSl) . we prove (b). 

Statement (c) about the Green matrix and the products of comparison solutions then follows 
from (IMl) and (1131). □ 



5 The diagonal of the Green matrix and discrete Agmon distance 

We next show that the approximations derived in ^for the solutions and Green matrix of by 
constructing a comparison (jl.2p are related to exact identities for Green matrices. In particular, we 
offer a discrete version of the discovery of Davies and Harrell in [TT], §4, that the diagonal elements 
Gnn of the Green matrix allow the full solution space and full Green matrix Gmn to be recovered 
formulaically. We build on significant earlier steps in this direction by Chernyavskaya and Shuster 
[71 [8]. As in jllj we furthermore point out connections between the diagonal of the Green matrix 
and an Agmon distance for p.ip . 

First observe that with respect to the diagonal of the Green function, there are critical differences 
between the continuous Schrodinger equations explored in [TT] and the discrete equations considered 
here. Consider that, due to Remark [H if 

Grnn ^min(m,n)^max(m,n) i^'^) 

is the Green matrix for some potential function Vn, then the same diagonal elements Gnn also belong 
to the Green matrix for an equation of type p.ip but with potential function Vn = —4 — V^- Thus 
the uniqueness of the representation of [11] is lost, at least to this extent. 

We have seen in ^ that for the comparison equation solved by the pair of functions (j4.ip , the 
factor Zn equals the square root of Gnn- Meanwhile, if Zn is given, then Sn is determined via (|4.4p . 
and consequently (j4.ip provides a basis for the solution space of the comparison (11. 2|) . and (14. 2p 
reproduces a full Green matrix for (|1.2p . Here we demonstrate that these implications do not rely 
on the construction of a comparison equation, but hold in generality for (jl.ip . 

Hence let Gmn be the Green matrix for any equation of the form (II. ip . and simply define 
Zn '■= VGnn- (If Gnn is negative, a canonical choice of phase could be assigned to Zn, but here we 
primarily consider the case where Gnn > 0.) We then use ()4.4p to define Sn for M + 1 < n < N 
viz., choosing the root analogously to ()4.6p . 



■= 2zz , ■ ^'-'^ 

Here we caution that this choice of the root of the quadratic equation for Sn will restrict the 
possible values of Vn in what follows. A pair of functions ip^ can now be defined by the ansatz 
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()4.ip . i.e., when expressed in terms of 

±1 

« 1 1 + Jl + Azlzl_A 
^t:=zn W V " • (5-3) 




Remarkably, with this definition, both Lp'^ and ip~ solve an equation of the form (jl.ip . where the 
potential function Vn is determined from Zn via 



n-l 



n 



1 + V1 + 444-1 



= ^^ii + ^ - 2, (5.4) 

provided that Vn > —2. (Else a different root must be chosen in (j5.2p .) To see that ip^ solve the 
same discrete Schrodinger equation, let us separately calculate 

^ = ^ + _ 2, (5.5) 

and note that since s|f' has been chosen to satisfy the equivalent of (|4.4p . the difference between 
these last two expressions is 

1 1 _^ 

4 4 

This leads to a theorem in the spirit of 

Theorem 5.1. Suppose that (II. ip /las too independent positive solutions for m < n < N , with 
N > M + 2, and denote the associated Green matrix Gmn- Since Gnn > for m < n < N , we may 
define Zn '■= \/Gnn- In terms of Zn, define Sn"^ and ip^ according to ()5.2p and ()5.3p . Then 

1. is an independent pair of solutions of (II. ip for m < n < N . 

2. Gnm = ZnZm HLm+l M < m< n < N . 

3. The potential function is determined from Gnn by a nonlinear difference equation, 

^ (Vl + 4G„nG„+ln+l + ^1 + 4GnnGn-ln-l) = (2 + Vn)Gnn- (5.6) 

Remark 5. The assumption that there are two positive solutions is related to the notion of dis- 
conjugacy in the theory of ordinary differential equations, cf. \17\ [1]. If, for example, > for 
n > Nq, then it is not difficult to show that no solution can change sign more than once, and that 
therefore the positivity assumption is satisfied for n sufficiently large. As will be seen in the proof, 
a necessary condition for the assumption is that Vn > —2. 

Per Remark [1] the positivity assumption can be replaced by the assumption that there are two 
solutions V'n such that {—l)'^Tp^ > 0. A sufficient condition for this is that Vn < —4 and a necessary 
condition is that Vn < —2. 
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Proof. The essential calculation was provided in the discussion before the statement of the theorem. 
Given that the Wronskian of ip~ and (/?"*" is 1, these two functions are linearly independent and 
therefore a basis for the solution space of 

(-A + yWV = o. 

Moreover, 

Gmn ''Pmin{rn,n)'^raa.-x.(rn,n) 

is a Green function for — A + Vn'^ . The crux of the proof is to show that Vn^ is the same as the 
original Vn of (jl.ip . 

\z\ 

Because Sn was defined such that 

- 4t = ^ . 

>— 'ri r 1 



we may rewrite (j5.4p as 



(5.7) 



From the definition of 2„ and the assumptions of the theorem, we know that for some independent 
set of positive solutions ipn of P-IP - with Wronskian 1, = ip^ipn- Therefore 



444±1 = 4('0n V'„±l)(V'nV'. 



2 



(V+V'^'il+V'nV'+il)'-!. 



Hence (j5.7|) yields 



2 + Ki^l = — ^ (^^C+l + V'n V'^+l + V'^^n-l + V'n ^n-l) 



2 + Vn, 



as claimed, and establishes (|5.6p . according to (|5.7p . □ 

Formula (j5.3p suggests that Sn can be related to an Agmon distance [21 [18], that is, a metric 
dAim^n) on the positive integer lattice such that every £^ solution of (II. ip satisfies a bound of 
the form 

and that as a consequence (;^~ decays rapidly as n — )• 00. Thus if Zn is bounded we expect an 

Agmon distance to be something like Yl'e=m+i '^i^'^ • write the Agmon distance in this way 
because a metric on the integer lattice must be in the form of a sum, as the triangle inequality is 
an equality.) In Agmon's theory, however, the distance function should be a quantity that can be 
calculated directly from the potential alone, and indeed, the estimates in ^ already imply some 
bounds of this form. As we shall now see, understanding the diagonal of the Green matrix allows 
the derivation of Agmonish bounds without the need to control expressions involving Vn+i — Vn, as 
in ^ We begin by showing that Gnn is comparable to {Vn + 2)~^ in a precise sense. 
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Lemma 5.2. Suppose that lim inf „_^oo ^« > C > and let Gmn be any Green matrix for (jl.ip . 
Define 

\/^+(c(C + 2))"+C(C + 2)- 

T/ien /or n sufficiently large, 



^{Vn + 2)(K_i + 2) + v/4 + 04 + 2){Vn-i + 2) ^ ^ v/(y^ + 2)(K-i + 2) + yM+TK + 2)(K_i + 2) 



(5.9) 

Remark 6. Note that the upper bound is of the same form as was found for the Liouville-Green 
approximation in (4.36). For a simpler bound Ka could be replaced in these inequalities by 



1 + -^^>Ka 
(see proof). 

Proof. The lower bound on Gnn is immediate from Statement (3) of Theorem 15. 1[ 



The upper bound in (15. Sp requires a spectral estimate. The Green matrix Gmn is the kernel of 
the resolvent operator of a self-adjoint realization of —A + V on £^([A^, oo)) for some A^, where the 
boundary condition at n = A^, + 1 is that satisfied by 99^. Since —A > on this space (as an 
operator), inf sp(— A + V)>C, and hence, by the spectral mapping theorem, ||(— A + F)~^||op < 
C""*^. Since Gnn = (cn, (— A + y)~^e„), where {e„} designate the standard unit vectors in it 
follows that Gnn < C~^. Inserting this into (j5.4p would already imply (j5.8p with Ka replaced by 
To improve the constant, replace only the terms Gn±in±i in (15. 4p by 1/C, getting 



1 I ^Gnn 

c 



(2 + K)<^^— . (5.10) 

^nn 

Since 

X 

is a decreasing function of x when x, y > 0, an upper bound on Gnn is the larger root of (I5.10p 
(which is effectively a quadratic). The claimed upper bound with the constant Ka results by keeping 
one factor T^i. + 2 in the solution of the quadratic, replacing the others by C + 2. 

The bounds on S'jf^ result from inserting the bounds on Gnn into (15. 2p and collecting terms. □ 

We can now state some Agmonish bounds. 

Corollary 5.3. Suppose that lim inf „_>.oo > C > and fix a positive integer m. Then the 
subdominant solution of (jl.ip satisfies 

(a) 



n 



J 00 
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(b) If, in addition, Vn^i — Vn £ then 



\e=m / 
Proof. Recall the representation (j5.3p . Because Zn is bounded, so is 



We then use the lower bound on S^^ from the lemma, but simplify by dropping the 4, which allows 
the product to telescope in a pleasing way, producing (a). 

For (b) we note that the additional assumption on Vn allows us to conclude that ip is well- 
approximated by the Liouville- Green expression in ^ Since z„ is again bounded, so is 



n 



using the ansatz (j4.ip . Finally, we recall (j4.6p . □ 
Thus when liminf„_>oo Vn > 0, a suitable Agmon distance dA{m,n) for (II. ip is given by 

n 

(HVi + 2) - In K a), 

or by 



^ + 2 + y/VejVe + 4) 

£=m+l ^ 

provided that £ The latter can be weakened to the simpler expression 

n 

l=m+l 

6 Some illustrative examples 

In Theorem 13.21 we consider the problem when sup„ < oo. Here we construct a potential V 
such that the boundedness condition of S„ is satisfied but S„ fluctuates as n — t- oo: 

Example 6.1 (bounded but fluctuating Sn)- Let V^ = V such that V [—4,0]. Then we may 
find a non-zero x G (—1, 1) such that 

x + - = {2 + V). (6.1) 

X 

The solutions to — A + V'^ = are given by 

(p- = x"" and (^+ = x-"". (6.2) 
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and the Wronskian W is x — x. Consider an asymptotically constant potential: 

= y + (-ifTyx2". (6.3) 

In other words, /3„ = /3n(j)n4'n — (— is summable and /3n{4't)'^ — Therefore, < 

suPn I n"=i(l ± /3i407)l < oo. For n > 1, 

n-l 

S„ = (-1)"+! + i-iyx^^) + (1 - (-l)"x2")S„„i (6.4) 

with 

Si = 1, T.2 = x'^-x^, S3 = l-x^ + x^-x^°. (6.5) 
Using ()6.4p . it is easy to prove that for k £ N, 

lim T,2k = and lim T,2k+i = 1- (6.6) 

fc— >oo fc— >oo 



Example 6.2. The main situation we have treated is where Vn — )• V^o ^ [—4, 0], with Vn — V^o G i'^, 
for which the solutions are of exponential type, with a subdominant solution. As a second case, let 
us suppose that — )■ 00 

Eigenfunctions that decay only polynomially are possible when Voo = or —4. Suppose, for 
example, that = n~" for some a > 0. This is a solution to a discrete Schrodinger equation with 
a potential satisfying 

Using a Taylor expansion, we find that 



k + lj \k-l 



y,= (^ = ^^(^ + 0(fc-). (6.8) 

Thus polynomial decay can be anticipated when the potential decreases like 'yk~'^. 
Corollary 6.1. Suppose that for some 7 > 0, 

where kW^ S ■ Then equation has a subdominant solution ipj^ such that 

lim /kKi+v'TO)^^ = 1. (6.9) 

For any solution ipk that is linearly independent of ip'j^ , 

^i(l-v^I+4^)^^ (6.10) 

converges to a finite, nonzero value. 

Next, we provide an example such that |</)+(/>~| is not bounded, yet the quantity J„ E and 
therefore Theorems 13.21 and 13.31 apply. 
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Example 6.3 (sparse perturbation). Consider a potential such that 

2 

n> 1 



:= Un + l)in-l) ' (6.11) 

l-i n = l. 

It is easy to verify that 0^ := 1/n is a solution to the equation —A + V^ = under the convention 
that (j)Zi = 0. By Corollarv 12.81 (p^ obeys |(/)+(/)~| ~ Cn for some constant C. 
Consider a potential V which is a sparse perturbation of V^: 

{W , 
V + ^ n = 2'' for some A: G N; ,^ ^ ^, 

(6.12) 
V otherwise. 

Under such definitions, /3n4>^(j)~ ~ C/n is sparsely distributed at powers of 2 and hence summable. 

Finally, we provide an example for which Liouville- Green approximation is accurate, while the 
potential fluctuates and diverges as n — )• oo. 

Example 6.4. Let be defined such that 

In" if n is odd; 
K=\, . «>2. (6.13) 

II if n IS oven, 

Then 

2Z . .1/2 ( ..3/2 + . .3/2 ) = W2~^ ^ \(n + 1)W2 (n - 1)^2 ) < (^-l^) 

n \ "^n+l n is odd n is even ^ ^ ' ^ 

Appendix: Second-order difference equations and orthogonal poly- 
nomials 

In this section, we will show how the discrete Schrodinger operator relates to orthogonal polynomials 
on the real line. We begin by recalling some standard facts; the reader may refer to |25^ [26] for a 
comprehensive introduction to the subject. 

Let be a non-trivial measure on M such that for all n S N, the moments are finite. In other 
words, 

\x\^d[x{x) < oo. (6.15) 

We form an inner product and a norm on L^(M, d^) as follows: for any f,g G L^(]R, d^u), we 
define an inner product and a norm as follows: 



if, 9) 



/(x)5(x)d/i(x), ll/f = / fixfdf^ix). (6.16) 

Jr 



By the Gram-Schmidt process, we can orthogonalize 1, x, x^, . . . and obtain the family of monic 
orthogonal polynomial on the real line with respect to the measure /i, which we denote as (P„(x))^q. 

For example, if ^ = V^tt then we obtain the Hermite polynomials; and if /i = X[-i,i]dx, 

then we obtain the Legendre polynomials. 
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Let {pn{x))'^=Q denote the family of normalized orthogonal polynomials, i.e., ||pn||^ = 1- It is 
well-known that the monic and the normalized orthonormal polynomials on the real line satisfy the 
following recurrence relations 

xPn{x) = Pn+l{x) + bn+1 Pn{x) + alPn-i{x), (6.17) 

xpnix) = a„+i(x)p„+i(x) + bn+iPn{x) + a„p„__i(x). (6.18) 

Note that ()6.18p above can be expressed as follows: 



(h 


ai 








Ol 



























( 1 \ 

pi{x) 

P2{x) 

■■■) \ ■: J 



( 1 \ 

pi{x) 

P2{x) 



(6.19) 



The tridiagonal matrix in (j6.19p above is called the Jacobi matrix. The recurrence relation (|6.18p 
can also be expressed in terms of the 2x2 transfer matrix ^n+i (x) as follows 



Pn+\{x) 

an+\Pn{x) 



On+1 



X - 



Pn{x) 
anPn~\{x) 



,n > 0. 



(6.20) 



l(^) 



Observe that the discrete Schrodinger operator with potential V and energy E on f can be 
written as 

-fn+l-fn-l + {Vn + 2)U = EU. (6.21) 

Compare (|6.2ip with (|6.18p . Note that the discrete Schrodinger equation (|6.2ip can be seen as 
having a„ = 1 and = + 2 and X = E. Hence, orthogonal polynomials associated with the 
measure with recurrence relations an = 1 and 6n+i = Vn + 2 evaluated at x = E can be seen as a 
solution of the difference equation ()6.2ip with initial condition {pQ{x),aop-i{x)) = (1,0). 

The solution to (j6.20p with initial condition (0,-1) (i.e., n = 0) are known as orthogonal 
polynomials of the second kind, {qnix))^=Q, where qn{x) is a polynomial of degree n — 1. Therefore, 
{pn{x))^^Q and {qn{x))^^Q form a basis for the solution space of the difference equation (j6.2ip . 

However, for the Schrodinger equation, we impose the condition that the solution is square 
summable (i.e. in £^(N)), a property that is not necessarily satisfied by pn{E). In fact, for any 




^l{xQ). 



(6.22) 



Hence, {pn{E))n is a solution if and only if £^ is a pure point /i. 

Recall the second-order difference equation (|1.13p studied by Geronimo-Smith [13] which was 
briefly discussed in Section [TJ Note that (|1.13p can be written in terms of a transfer matrix 



y{n + 1) 
y{n) 



d{n+l)' 



q{n) 
d{n+l) 







y{n) 
y{n - 1] 



(6.23) 



which resembles the transfer matrix An+i{x) in (|6.20p . Hence, techniques developed to study the 
asymptotic behavior of orthogonal polynomials can be applied to study ratio asymptotics of the 
solutions, which determines whether the limit lim„_5.oo y{T^ + l)/y(^) exists and what the limit is in 
the case that it does. For (|1.13p and given that y{n) = nj=no ^0')' I'^-tio asymptotics means 



y(n + 1) _ u{n + 1) 



y{n) 



u{n) 



(6.24) 
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which explains why it was reasonable for Geronimo-Smith to assume that hm„_^oo^(^ + 
exists should the convergence rates of q{n) and d{n) be sufficiently fast. 

For the asymptotic analysis of Pn{x) by means of the transfer matrix when the coefficients are 
asymptotically identical (meaning an ^ a, bn ^ b), the reader may refer to [32] . 
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